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Abstract 

The problems of neutrino production states, prepared to the oscil- 
lation process, in the case of non standard interactions, are briefly dis- 
cussed. Quantum neutrino states are determined from the dynamics of a 
production process. We show, that even in models where only left-handed 
neutrinos are introduced, the standard adopted procedure is valid only ap- 
proximately. Entanglement between neutrino masses or between masses, 
flavour and spins cause, that their quantum states are mixed. 

1 Introduction 

Without the Non-Standard Interaction (NSI) the existing theory of neutrino 
oscillation works very well. In the Standard Model (SM) neutrinos interact 
through the left-handed vector currents. Then the relativistic neutrinos (an- 
tineutrinos) have only negative (positive) helicity. The Maki-Nakagawa-Sakata- 
Pontecorvo (MNSP) mixing matrix present in the model and definite helicity 
cause that neutrinos are produced (and detected) in a pure Quantum Mechani- 
cal (QM) states. The neutrino helicity is only a "spectator", it does not change 
at any phase of an oscillation process, even if a flavour is changing. 

\v a ,D = £t£><4>, KA) = J>«,<fa,t>. (1) 

i i 

The fact that the NSI can change the process of neutrino oscillation for the 
first time was observed in pQ, and then has been considered in various physical 
processes (see e.g. [2j |3l H] and references therein). In all this approaches 
instead of the states (TjQ) the new one, connected in some way with dynamics of 
the processes are defined separately for production (s) and detection (d) 

K> = l^> = £^>*>' (2) 

i i 

where the mixing matrices i and Ui i depends on the NSI 

\U^\ 2 oc\(^;f P \H p \l a ;i P )\ 2 , lUlf oc \(l^f D \H D \ui;i D )\ 2 . (3) 
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with H p (H D ) being the production (detection) Hamiltonian. 

Such description of the oscillation phenomena in frame of NSI is widely ap- 
plicable (see e.g. [SJ|I)1[7]). However, it was observed before (see Ref. 8, 9 ), that 
such approach is valid only approximately. There are models which require an 
extension of the traditional approach presented above. There can be found many 
examples of such theories in literature, where such an extension is necessary. Let 
us mention only a few of them: models with two Higgs doublets [TO], Zee-Babu 
mo del [IT], models with extended symmetry group[jT2" | ll3 | or supersvmmetrv[T4"] . 

In this work we would like to reconsider once more the impact of NSI for 
production process. We concentrate on processes which are interesting from the 
practical point of view as for example; pion decay used in the present accelerator 
experiments, nuclei beta decay for planed future beta beams, and muon decay 
for future neutrino factories. 

We also discuss the case, where only left-handed neutrino fields enter a 
model and relativistic neutrinos (antineutrinos) always have negative (positive) 
hclicity. Then the difference between conventional approach (pure neutrino 
state) and the real one (mixed neutrino state) are not connected with obvious 
mixing of states caused by the different neutrinos helicities . 

As we will see in the next Chapter, the mixed state, if appear, is connected 
with subtle things of an entanglement between neutrinos or antineutrino in 
different mass states (in the muon decay case, where neutrino- antineutrino ap- 
pears) or between masses, flavour and spins, when only one neutrino is produced 
(pion decay, beta decay). 

2 Quantum state of neutrinos produced in the 
processes with NSI 

The decay of pions, nuclei, or muons are in the lowest order described by the 
general d=6 effective Lagrangian 

£ NS i = -2v / 2Gf £ tr,@aT A Psi>b)@ c TAP v i>d) +h.c, (4) 

A = S,V,T 
e,ri=L l R 

where the field tpaijpb) describes the produced neutrino (antineutrino). In many 
extensions of the Standard Model (MS) only left-handed neutrino fields ip a ,L 
are present, in such a case the Lagrangian (Eq.Q) simplifies (five parameters; 
£ R -R(L)' £ fi r(l) an d e lj l are then not vanishing). If two left-handed neutrinos 
appear (e.q. in addition d = v 1 like in the muon decay) then only one scalar 
and one vector term remains (e R L , e\ L ^ 0). Then the Lagrangian (UJ), by the 
Fierz rearrangement theorem (see e.q.|15)V has the form similar to the one used 
for studding the NSI in matter (see e.g. [T51 ITT1 ITS] ). 

£ NS i = -2y/2G F £ % C ^PcTfPm) (jluPcf) +h.c, (5) 

C = L,R 

where e\ L =>■ -e^'Jj L , e s R L =>■ 2e£jg R and i\) a =>■ v a , ip b /' , ipc => /, ^d^v p . 
On what follows we will consider this form of the NSI. 

Despite that the main effects of this interaction is to modify the neutrino 
oscillations in matter, it is also claimed that NSI modifies the production and 
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detection states such that those states are given by (not normalised) see e.g. 

EEgngi: 

K) = K) + E ^l^)' ( 6 ) 

(>4\ = E £ U^\, (7) 

a— e,/i,r 

where the parameters e s ^ are connected with the complex coefficients e^Jg C in 
the way which depend on the given production (s) and detection (d) process. 

We would like to reconsider once more the impact of NSI given by (0) for 
the production process. In order to calculate the state of neutrino produced 
in some reaction we write the density matrix describing a final state of this 
reaction [SJ H] and then, using this final quantum state, we take the partial 
trace over unobserved particles in order to obtain the state of neutrino which 
we are interested in. We will also work in a mass base rather than in flavour 
base. The mass base is convenient because it is well defined and suitable for 
considerations of the oscillation phenomenon. To move from flavour to mass 
base we assume as usual that flavour neutrino is given as a linear combination 
of mass eigenstates i.e. 



V,x = E UaiVi (8) 
i 

where U is some unitary matrix, in SM it is just the MNSP mixing matrix. 
Generally left- and right-handed neutrino fields can have, depending on the 
neutrino mass model, different mixing matrices, but in our case Eq. © is 
sufficient. 

The modification of initial and final state given by eq. ([6]) and assumes 
that the all spin amplitudes are added coherently. This is not true in general for 
NSI ||SJ). Therefore we will consider both cases i.e. coherent and incoherent NSI 
contribution to the production process. Furthermore the situation gets more 
complicated with two neutrinos in final state, therefore we will consider a muon 
decay as an example of production process 

M~ -> eTviVj, (9) 

but to simplify the calculation we assume that neutrinos are Dirac particles. 
Now in © the "i" and "j" indices indicate the neutrino mass. 

If neutrinos are left-handed the most general effective Lagrangian (with the 
SM part included) in the mass base for the process has the form 

£ M = -2V2G F [gfj {ViP R e) (JIPlUj) + g% faj a P L e) (Jl^PlVj)] + h.c. (10) 

Taking into account the SM part we should replace the e's factors defined by 
eq. © by 

£ Tf] L G "/3 = + S^p , (H) 

then the relation between Lagrangians (0 and (|10p are given by 

al = E Ki^Upi - u* ei u M , (12) 
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and 



gf j = 2Y f Ki<fu Pj . (13) 

We will consider two cases. First we assume that in (J5j) only one term C 
= L appears, so 

£ af ^ and £ ^ R = °- In tllis 

case the NSI contribute 

to the same helicity amplitude as VF ± boson exchange in SM (coherent case). 
The density matrix which describe the antineutrino - neutrino state produced 
in muon decay © in the mass-helicity base (it does not depend on an initial 
muons polarisation, and whether the final electrons polarisation is measured or 
not) has then the form 

Q(i,j;k,l) = x*jXk,h (14) 

g v I 

where Xi,j = an d N — j I 9ij I 2 - F° r tne left-handed and to a good ap- 

proximation massless neutrinos, their helicities are unambiguously determined, 
and therefore in formula ([T4"]) their are not indicated explicitly. Such density ma- 
trix describes pure QM two neutrinos state (Tr(g 2 ) = 1), which can be written 
in the form 

v») = E(^-)i^) = jfiM +£(C/n^»- (is) 

id a,P 

In order to get one particle state, we must take the partial trace over the second 
particle degrees of freedom (g"(i;k) — J2j Q(h3'>k,j)) for antineutrinos and over 
the first mass indices (g u (j; I) = Q(i,j', h 0) f° r neutrinos and we obtain 



g ^ k >~ rr(GGft) ' flJ|l) " Tr(GW) ' 

where U and G matrices arc defined in Eqs © and (jlip respectively. Generally 
such states are not pure, because 

, Tr[&G&G] 
Tr{Q } = (Tr[Gt G ])2 * 

Only for very specific NSI couplings such that Tr{A 2 ) = (TrA) 2 with A = G^G 

both neutrino states produced in muon decay are pure. In general these state 

can not be written as a pure state because the NSI caused the entanglement of 

neutrino and antineutrino produced in muon decay. Only if the initial state (fT5")) 

is factorisable then the produced neutrino can be described as a pure quantum 

mechanical state. If we look at the models beyond the SM, then the parameters 
f f'c 

e a g are usually parametrised in the way 

r ff'C _ h »f h h' h o\ 

where h a f and A are some couplings and a scale of the New Physics. Unfortu- 
nately this not guarantee that the state (TT51) will be pure. This will be the case 
if h ae = d ae or = 5 all . In this situation the state © will be the correct 
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one. In general the difference between © and (fTTT) is only of order e 2 so then 
with the present bounds on NSI [T7J [T5] the error made by taking a pure state 
instead of the density matrix is very small and the standard approach is justi- 
fied, but it must be remembered that © is only the first order approximation. 
This situation is typical in reactions which involve more than one neutrino in 
the final state and the NSI giving a coherent contribution. 

Now we consider a general case. Let us assume that NSI contributes to the 
scalar and vector part in (JTUJ) so e e ^ 7^ and 7^ 0. Such interactions 

can, for example, appear in theories with charged scalar boson like models with 
more than one Higgs doublets [20]. The vector and scalar terms from (fT0|) 
contribute to the same helicity muon decay amplitude (if m e 7^ once more the 
NSI contributes coherently) but with different kinematical factors giving finally 
the density matrix for muon neutrino 

q = ^(Bf(E) + Cg(E)~2Re[(g v ) T (g s y}(^)h(E)), (19) 

where f(E) = 6M - 8E + O(^), g(E) = 3(M - 2E) + O(^gf-) and h{E) are 
functions depending on neutrino energy E in muon rest frame and muon mass 
M. The matrices A and B are defined in the way B = (g v ) T (g v )* and C — 
(9 ) T (9 )*■ Such density matrix is linear in NSI parameters (third term in (|19|l ). 
but unfortunately this term is proportional to small factor (t^) and therefore 
we neglect it. Then the normalising factor is given by 

N = J(E)TrB + g{E)TrC. (20) 

f f' L f f 1 R 

When neglecting the electron mass e A and e a/3 contribute to the different 
hclicitics amplitudes (incoherent case). State given by (fT§|) is mixed, Tr(g 2 ) is 
given by 

Tr(g 2 ) = -^(f 2 (E)Tr(B 2 ) + g 2 (E)Tr(C 2 ) + 2f(E)g(E)Tr(BC), (21) 
N 

To get neutrino in a pure state for any value of neutrino energy, couplings 
should satisfy very restrictive conditions Tr(B 2 ) = (TrB) 2 , Tr{C 2 ) = (TrC) 2 
and Tr(BC) = Tr(B)Tr{C), which do not follows from any symmetry of the 
considered models. 

We see that by taking a pure state instead of the density matrix, as an initial 
neutrino state we overestimate the impact of NSI because it is of order e 2 at 
most and there is no way to obtain the state (]6]) as some kind of approximation 
of formula p^|) . In (IT9")) there is also term linear in e but it is proportional to 
the small ratio ^ and we neglect it. 

What can we say about the quantum states, when only one neutrino is 
produced, as in the pions decay or beta decay of nuclei? As we discussed 
previously - as many as five operators describe production of the left-handed 
neutrinos, but there is no mass entanglement between them. As as example let 
us take effective Lagrangian for antineutron beta decay without tensor operator 



g? L ' R {viP R e) {nP L , R p) + gY L ' R (vil a P L e) (n<y a P L ,R p) + h.c. (22) 
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If the coupling have similar structure as before (e.g. gf L — J2 a 9a L ^ai) then 
the state of produced neutrino is once more mixed which is caused by the decay 
amplitudes 

/(An, A p , A e , vi) = Ki9« /a (A„, A p , A e ). (23) 

a, A 

Where /A(A n , A p , A e ) are the amplitudes calculated for a given operator A = 
Sl,r,Vl,r- The state will be pure only if all /a's for any spin configuration 
of accompanying particle are equal (it can not happened in practice) or the 
couplings do not depend on the operators A (strong constrain). Once again, 
we have kind of entanglement, in this case, between neutrino masses, flavor 
structure and spin projections. 

We considered the NSI with left - handed Dirac fields only, as the possible 
values of the right-handed couplings are strongly constrained, especially with 
the requirement of gauge symmetry (see e.g |21j). In models which accept 
right-handed chiral fields, neutrinos (antincutrinos) with a positive (negative) 
hclicities are produced and their states are obviously mixed. 

3 Conclusions 

Quantum state of neutrinos, which are produced in the present and future os- 
cillation experiments was re-discussed. It was shown that these states depend 
in a crucial way on the details of the production mechanism beyond the Stan- 
dard Model. Assumption that they are always pure quantum states is incorrect. 
Even when the helicity of neutrinos is clearly defined, their mass states and 
flavours remain entangled, causing that the full state of a single neutrino is 
mixed. Sometimes the effect can depend linearly on the parameters of New 
Physics and when one is searching for such an effects in neutrino oscillations 
must keep this in mind. 
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